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Abstract 

m 

We construct a class of scalar field models coupled to matter that lead to the dependence 
\ of masses and coupling constants on the ambient matter density. Such models predict a 
deviation of couplings measured on the Earth from values determined in low-density astro- 
physical environments, but do not necessarily require the evolution of coupling constants 



with the redshift in the recent cosmo logical past. Additional laboratory and astrophysical 
tests of Act and A(m p /m e ) as functions of the ambient matter density are warranted. 
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1 Introduction 



Perhaps the most astonishing fundamental observation of the last decade was the discovery 
of dark energy. So far, all cosmological data are consistent with the simplest possibility: 
dark energy is just a new fundamental constant of nature, which does not evolve over 
cosmological redshifts. On the other hand, it is intriguing to think about alternative expla- 
nations associated with this profound change in infrared physics. The most straightforward 
way of implementing such a change is the introduction of a new ultra-light scalar degree of 
freedom associated with quintessence [1]. The coupling of this scalar field to matter may 
be the source of new cosmological phenomena such as an apparent "breakdown" of Lorentz 
invariance connected to the CMB frame, the existence of a "fifth force" mediated by scalar 
exchange, or a change of couplings and masses with time. Thus, the search for these exotic 
effects acquires a new actuality in the "dark energy" age. 

A potential hint on the difference between the laboratory values of the fine-structure 
constant a [2, 3] and the one derived from quasar absorption spectra at high redshifts [4], 
Aa/a ~ —0.6 x 10~ 5 (z ~ 0.5 — 3), has triggered a series of new phenomenological and 
theoretical studies of "changing couplings" [5]-[16]. Subsequent studies that employed the 
same "many multiplet method" have shown Aa consistent with zero with the same accuracy 
[17] (however, a critical discussion of this result is found in [18]). In a recent development, 
the comparison of H 2 spectra obtained in the laboratory and at high redshifts, yielded the 
first hint of a possible change of m p /m e at the three sigma level [19], whereas the analysis of 
the inversion spectrum of ammonium at z — 0.68 led only to an upper limit on A(m p /m e ) 
[20] . Null results for time evolution of couplings are reinforced by the chemical composition 
of Oklo rocks [21] and meteoritic abundances of rhenium, that provide stringent constraints 
and go back to z ~ 0.2 — 0.4 [22, 23]. Finally, a look into the much deeper past, to the 
time of the Big Bang Nucleosynthesis (BBN) [24, 25, 26], leads to the conclusion that the 
couplings could be different from present day values by no more than a few percent at that 
time. 

Despite the controversial status of the non-zero claim for Aa/a there have been a 
significant number of attempts [5]- [16], [27] - [29], to build simple models that could account 
for possible effects of order O(10~ 5 ). Theoretical models of the variation of other couplings 
were discussed in Refs. [30, 25]. The simplest Lagrangian that enables the variation of 
a was written down by Bekenstein [27], and involves the coupling of an ultra-light scalar 
field to the kinetic term in the electromagnetic Lagrangian, Qp^-F^F^ ', where M* is 
a large energy scale, comparable to the Planck scale and (p is a constant. An attractive 
consequence of this interaction is the absence of the evolution of the coupling in the radiation 
dominated epoch after the annihilation of electron-positron pairs [27] . The drift of a in time 
comes about later during the matter dominated epoch due to a nonzero, O(10~ 4 — 10~ 3 ), 
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electromagnetic contribution to the nucleon mass, which supplies a coupling of tp to the 
baryon energy density. The time evolution of a, including the sign, is calculable in terms of 
one free parameter, ((p/M*) 2 , and the result is a linear dependence of a on the logarithm 
of the redshift in the matter dominated epoch. Unfortunately, between redshifts z ~ 1 
and z — 0, Aa/a is predicted to be far too small in the minimal model of Bekenstein; 
several orders of magnitude below the modern observational capabilities, once the fifth 
force constraints on £p/M* are taken into account. In attempt to resuscitate this model, 
it has been suggested that a coupling to dark matter [28, 8, 9] and/or a self-interaction 
potential [10, 11] of the scalar field can drive its evolution in a much more efficient way 
than the baryon energy density. Inescapably, these models have more parameters {i.e. a 
nearly arbitrary potential V(<p)), which results in a loss of predictivity. Thus almost any 
redshift dependence for a(z), including oscillatory behavior, is possible, and no concrete 
predictions can be obtained from these theoretical models. If the experimental limits on 
Aa are enforced, it has proven to be difficult to remain consistent with a substantial change 
in the fine structure constant at z ~ 1 and satisfy Oklo and meteoritic constraints. Still, 
fairly large classes of models (see e.g. Refs. [12], [14]- [16]) are known to pass these 
requirements. On the downside, if V(ip) supplies the driving force for the scalar field, there 
are no convincing arguments why ip would not evolve deep into the radiation domination 
epoch resulting in vastly different couplings during the BBN. In addition to these specific 
difficulties, there is a generic naturalness problem that all such models suffer from: it is 
difficult to argue that the mass of tp, as well as other terms in V(ip), are protected to 
the scale of 10 -33 eV to allow for the requisite late evolution, and at the same time allow 
the 93-field to have significant couplings to matter fields. Any conceivable loop effect or 
nonperturbative QCD vacuum condensate would tend to induce mass parameters about 
twenty orders of magnitude in excess of 10~ 33 eV (See e.g. [31]). 

In all models of changing couplings discussed so far in the literature, the temporal change 
of a dominates over the possible spatial variation of a. Although some suggestions were 
made that the spatial variation may dominate over temporal variations [32], no explicit 
models where found to date. In Bekenstein-type models, the spatial variation of couplings 
is caused by matter inhomogeneities and thus follows the profile of the gravitational field. 
If the coupling of the scalar field to matter is not stronger than the matter-gravity coupling, 
one expects the change in <p to be less than the variation of the metric. For example, the 
difference between coupling constants on the surface of the Earth an in orbit is not going to 
exceed 10 -10 and in practice will be much smaller once fifth force constraints are imposed. 
In this paper, we show that such a conclusion is not generic, and there is a whole class of 
models where spatial variations can be more pronounced than the cosmological variations 
in the recent past, opening new possibilities for searching for Aa and Am as functions of 
the matter density. 
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A key to this proposal is to choose the couplings of a scalar field to matter to be much 
stronger than gravitational. At first glance this would appear to only worsen the problem of 
a fifth force induced through scalar exchange. However, this might not happen if the matter 
density itself leads to the effective suppression of the linear scalar field coupling to matter 
[33, 34], and/or the range of the scalar-mediated force becomes shorter than the one needed 
for conventional fifth-force experiments [35, 36, 37, 38]. Models that escape prohibitive fifth- 
force constraints may predict spatial variation of a and m p /m e that exceed recent temporal 
variations. As we will show in the remainder of this paper, such constructions can be 
achieved if the dynamics of the scalar field in low-density environments is determined by 
its self-potential, while in regions of large overdensities, the dynamics of ip is set by its 
coupling to matter. In a large subclass of such models where the matter-^ coupling is 
much stronger than gravitational, the temporal evolution happens on the time scales that 
are much shorter than cosmological. In these models, the global temporal evolution of the 
scalar field could be finished a long time ago, and on average the "cosmological" values of 
masses and coupling constants remain constant in space and time. We denote the values of 
couplings and masses in rarefied low density regions as {a r , m r }. The value of the coupling 
constants on Earth or in any high-density environment can be different and are labeled as 
{ad, rrid}. In the next sections we discuss the essential features of the model, determine 
observational constraints on its parameters and argue that independent tests of ad — a r 
between different points in space are indeed warranted. 



which can be viewed as a generalization of a scalar-tensor theory of gravity. In this expres- 
sion, Mpi = (87rGAr) _1 / 2 is the reduced Planck mass, is a dimensionless scalar field with 
M* being the analogue of the Planck mass in the scalar sector. The functions Bpi(4>) give 
the ^-dependence to the gauge couplings in Standard Model (SM), and the sum is extended 
over all SM gauge groups, ipj represents Standard Model fermions that are coupled to <j) 
via the functions Bj(<f>). After performing a ^-dependent rescaling of the matter fields, one 
is allowed to remove the ^-dependence of the kinetic terms for the SM fermions ipi and 
keep only couplings to the mass terms. If needed, the interaction (2.1) can be generalized 
further to include Higgs bosons, cold dark matter particles, etc. 



2 Scalar field models of a(p) and m(p) 



The starting point for our analysis is the matter-gravity-scalar field action, 




(2.1) 
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Among couplings to the SM model fields, the couplings to quarks, gluons, photons 
and electrons are the most important. At lower energies, we can abandon the quark-gluon 
description in favor of an effective coupling to nucleons and reduce (2.1) to a more tractable 
form, 

r ( M 2 M 2 

_Bf^_ f ^ u+ £ /i ; (o)/,/,r,r, ). (2.2) 



j=n,p,e ) 

Since we are going to consider couplings of that are essentially much stronger than 
gravitational, the stability of the model will require that V(<f)) and the -B«(0) functions 
have a minimum with respect to <p. In what follows, we shall adopt the following ansatz, 

V(cf>) = A + ^A 2 (0 - O ) 2 + 5,(0) = 1 + \ti{<t> - 0*) 2 + -, (2-3) 

where ellipses stand for cubic, quartic etc. contributions around the minima. Here £j, O 
and (pi are arbitrary dimensionless numbers; A and A 2 have dimensions of [Energy] 4 and 
we are tempted to choose A to be equal to the current dark energy density to "solve" the 
dark energy problem. 

A further simplification of the quadratic ansatz comes from the assumption that the 
proton and neutron -E> p ( n ) functions are mostly induced by the gluon 5-function, and thus 
are approximately equal. With these simplifying assumptions, we can take 

K-4> V = 0m] £n ^ £p = 1; 00 = 0. (2.4) 

The normalization of £ p ( n ) to one can be attained by rescaling M*. In principle, a negative 
value for £ is also possible, but in this section we shall restrict our discussion to positive £'s. 
Of course, the relations (2.4) are only approximate, and possible violations at the ~ 1 — 10 
per mill level are naturally expected due to the nonzero quark and electromagnetic content 
of nucleons. The choice of O = can always be achieved by a constant shift of 0. The 
ansatz (2.3) and (2.4) is very similar to the D amour- Polyakov model [34] (see also [33]), 
where all couplings to matter fields exhibit the same minimum. In the same vein, we assume 
the same minimum <p m for Bp(4>) function. There are two important differences in our 
approach compared to the Damour-Polyakov models: we take M* to be much smaller than 
the Planck mass, and introduce a self-interaction potential that has a different minimum 
than the minimum of the -B«(0) functions. 

In this section we disregard higher-order nonlinear corrections to V c g, postponing their 
discussion to Section 4. Furthermore, we assume a region of relatively uniform matter 
density p. In such regions, the scalar field equation of motion takes the following form 

M, 2 D0+^ = O, (2.5) 
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where the effective potential is given by 



Kff = A o + ^A 2 2 + ^(0-0 m )V (2.6) 
This potential creates the minimum for the scalar field at 

0min = 4>m ~T~ i ( 2 -7) 

P + A 2 

and the physical (canonically normalized) excitation tp around this minimum has a mass 



By definition, the longest range for the (^-mediated force is achieved in vacuum at p = 0. 
It is instructive to present a numerical formula for \ c s at p ^> A 2 : 

, M* /10 24 GeV cm~ 3 \ 1/2 
A cfr = 7 x ID" 3 cm x ^ j , (2.9) 

which shows that for an extreme case with a weak-scale M* and terrestrial matter densities 
the range of the force falls under one millimeter. 

If the spatial extent of the mass distribution is much larger than the Compton wave- 
length of the physical excitations of 4> the effective interaction with a "test" nucleon takes 
the following form, 

r AT AT ( \ I <^ A 2 ^ V 2 \ , ot M 

from where we can read a p-dependent mass of a nucleon, 

m - = m " i 1 + 2(frt) • (2 - n) 

and the scalar-field-corrected Newtonian interaction potential between two nucleons sepa- 
rated by distance r, 

V(T) = (l + exp(- w ) x 2 -§£^_) ■ (2.12) 

Perhaps the most interesting case to consider is A 2 3> p for low density environments, 
such as e.g. the interstellar medium, and A2 C p for high density environments such as 
stars and planets. In that case, the change in the nucleon mass and the fine structure 
constant can be expressed as 

Am at m Nr - m Nd 2 



m . 



m N m N 2 

Aa a r - a d £ F 4> 2 m 



(2-13) 

a a I 

and we assume that 4> 2 m and t^F^m are m uch less than one. Notice that £ F can be as large 
as £ F ~ 0(100) without violating the assumption that £ n ~ £ p ~ 1. 



5 



3 Experimental constraints on the model 



All experimental constraints on the model described by (2.10) can be divided into two broad 
categories. The constraints coming directly from the quadratic couplings of to matter 
to a large extent do not depend on the position of the minimum of and on whether this 
minimum can be reached for a realistic size of an overdensity in question. The second group 
of constraints follows from the linear coupling of to matter, which are very sensitive to 
the position of and on the size of the overdensity. 

General remarks on chameleon-type models. Before we proceed with the analysis of our 
model, we would like to make several remarks regarding chameleon models. It has been 
shown that an appropriate choice of self-interaction potentials [35] relaxes gravitational 
and astrophysical constraints on the density dependent interactions of the scalar field [39] 
caused by a shift in the field 0. Density dependent couplings could in principle cause a 
shift in the field value as well as its mass. However, the magnitude of the shift will depend 
sensitively on the local density and the length scale over which the shift occurs. For the 
chameleon mechanism to work, we must require that dV e g/d<j) S> M 2 R~ 2 5<p or m 2 R 2 ^> 1, 
where R is the characteristic scale of the density enhancement and <50 is shift in from 
the low density (cosmological) solution of to the local value. When this condition is 
satisfied, the field and its mass inside the overdensity and away from the boundary will be 
determined by the solution of a spacially homogeneous equation dV e ^/d(j) = 0. 

To begin with, let us consider models with a linear coupling to density 

V eB = V(<p)+0fxp, (3.1) 

where we use the canonical normalization for the kinetic term of <p and choose f3 = M~ l . 
In chameleon models with a quartic potential, V(ip) = A<p 4 , the shift in <p is tiny as long 
as A is not tremendously small. The effective mass when dV eS /d(p = is A 1/6 (/3p) 1/3 and 
for A ~ 0(1) is far greater than R^ 2 - In this for all stiff potentials, the local en- 

vironment determines the field dependent couplings. Another example is the quintessence- 
like potential V((p) = M 5 /(p with M ~ 10~ 24 M P i. The cosmological background solu- 
tion for p ~ Mp/ gives (p ~ M P1 so that the potential and cosmological density are of 
the same order of magnitude. Naively, the local solution of dV e a/dip = would yield 
<p 2 = M 5 M P1 /p e ~ 10 _27 M|! taking p e ~ 10 _93 M^. However in this case, the gradient 
term in eq. (2.5) dominates and the correct solution for the local value of is a small shift 
from the background value of order 5<p ~ p® R 2 e /M Pl or S(p ~ gR®M Pl ~ KT 9 M P1 where 
g is the local acceleration on the Earth's surface. In this case, the chameleon mechanism 
is not operative. It can be restored if one adds a constant, M 4 to V(ip). By doing so, the 
mass scale, M, can be made significantly smaller M ~ 1CT 30 M P1 [36] and now the back- 
ground solution for <p yields a very small value <p ~ 1CT 15 M P i. In this case, the gradient 
terms can be safely neglected and the local solution for <p is indeed given dV e s/d(p = 0. A 



6 



similar argument can be made to show that the local density has virtually no effect on the 
background field value for potentials of the form V(ip) = — /i 4 ln(<p/M P i) with p ~ 10~ 3 eV 
that was discussed recently in [40]. 

Exact solutions for spherical (under) overdensities . The model based on the effective 
potential (2.6) considered in this paper is simple to analyze, as the field equations are 
linear. This allows us find an analytic solution for a spherical region of constant density, 
and consider both large and small M* limits. The general form of the solution for as the 
function of radius for a spherical region of density p\ of maximal extent R surrounded by 
the infinite region of density p 2 takes the following form: 



0i + — sinh(mir), r < R; 
0(r) = { ' (3.2) 



^2 + 7 exp(— m 2 r), r > R, 



where the constants of integration 

A = (4>2-4>i)(l + flm 2 ) (3 3) 

m 2 sinh(mi-R) + m x cosl^mxi?) ' 

B = exp(m 2 i?) x ^ - M^^R) - Rm, cosh(rr^)) 

rri2 sinh(mi_R) + mi cosh(miit) 

are uniquely determined by boundary conditions, for which the physical choice is 0(oo) = 2 
and 0'(O) = 0. In these formulae, mi, m 2 , and <f>i, <p 2 are the mass and vacuum expectation 
values calculated according to (2.8) and (2.7) for p = p x , p 2 . Eq. (3.2) is a generalization 
of the solution previously found in Ref . [39] . 

Adapting this solution to the case of p\ = p e and R = R 9 , it is easy to see that in the 
limit of large M*, or more precisely miR^ ~ p^ 2 'M~ l i? e 1, the solution simplifies to 
a quadratically rising function on the inside and a 1/ r-falling function on the outside (for 
simplicity, we also take A 2 , p 2 — > leading to 2 = 0), 



0(r) = < 



f -0i - ^f) , r< R; 

(3.4) 



Pi-t-i r>R. 



It is easy to see that this solution exactly follows the gravitational potential profile, with 
0i = 0m- The change in the nucleon mass between r = i? e and spatial infinity is 

— — = -(0 (oo) - (i?®)) ~ {gR e ) x (10 - 10 ) x M2 . (3.5) 

which is exactly what one anticipates in a linearized scalar-tensor theory of gravity with 
the relative strength of spin-0 to spin-2 exchange given by 2M| 1 0^ l /M^. Furthermore, this 



7 



limit, namely p e M~ x R e <C 1 is guaranteed in our model as long as M* > 10 13 GeV. This 
regime is of no further interest to us in this paper, as it has been investigated in a number 
of previous publications. 

The opposite regime is achieved when rriiR ^> 1, in which case the interior solution 
quickly adjusts to = 0i for r < R. Moreover, if in addition m2-R ^> 1, the spherical 
symmetry of the problem becomes irrelevant, and the solution degenerates into being 
frozen to its respective minima, <pi and 02, everywhere in space except for a small region 
near the surface separating two density regions. Defining this surface as z = 0, we can 
write down a simplified form of (3.2) which now takes the form: 

<f>{x) = (3.6) 
I 02 + %^exp(-m 2 z), z>0. 

A smooth boundary would lead to an adiabatic adjustment of between its minima provided 
that the Compton length of the scalar is much less than the characteristic scale of density 
change. For the atmosphere this scale is given by A atm = l/\dlog(p a tm/ P®)/dz\ ~ 1km, 
and the adjustment of the scalar field will occur as long as 

Acflf (0aim) < A otm , (3.7) 

ensuring that in a terrestrial laboratory environment is exponentially close to m j n (2.7). 
Neglecting the small A2-proportional contribution to the scalar field mass, we find that the 
condition (3.7) implies that 

M* < 1km x pill M* < 10 9 GeV (3.8) 

Clearly, this choice of M* selects scalar models that are significantly more strongly coupled 
than Mpi-normalized models and therefore the behavior of on Earth will depart drastically 
from the gravitational potential ~ goo ~~ 1- Next we explore whether the choice of strong 
coupling (3.8) can survive gravitational and astrophysical constraints. 

Astrophysical constraints. First we discuss the astrophysical constraints on the model 
which employs a 2 coupling to photons and nucleons (the linear coupling is suppressed by 
(A 2 /pi) 2 and is assumed to be 1). It is clear that the quadratic coupling will be less 
severely constrained than a linear coupling by the thermal emission rate of 0-quanta from 
the hot interiors of stars. Indeed, the overall emission rate scales as M~ 4 rather than f~ 2 
as one would routinely find in an axion-type model. As a result, instead of a lower limit to 
f a or order 10 9 — 10 10 GeV, we expect to find a much more relaxed bound on M*, of the 
order of the electroweak scale. 

Let us calculate the emissivity of quanta due to pair annihilation of photons. The 
amplitude for this process is induced by the t^p^F^F^ term in the Lagrangian, leading 
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to a cross section for this process in the center of mass frame, 



£ 2 ^ 2 

This cross section results in an energy loss (Energy/volume/time) for a thermalized gas of 
photons at the level of 

r 77 _^ = n^ua^) = ^- ~ 0.06 x (3.10) 

Comparing this to the typical limit on T = e x p core < 10~ 14 MeV 5 that follows from the 
constraints on the emissivity of light particles in cores of supernovae [41], 

< 10 19 erg g" 1 s- 1 at p core = 3 x 10 14 g cm" 3 , T core = 30MeV, (3.11) 

we obtain a typical sensitivity to the coupling of to photons, 

M^p 1 ' 2 Z 3TeV. (3.12) 

This limit is admittedly not very precise, as it is quite sensitive to the temperature of the 
core, and more conservative assumptions about T core may result in relaxation of (3.12) by 
a factor of a few. Other channels of 0-production from light species, such as — > e00 or 
e + e~ — > 00 will be further suppressed by the smallness of electromagnetic couplings or by 
the ratio of m e /T. 

The constraint (3.12) is not far the limits on the 6-dimensional Planck scale M 6 in models 
with two large extra dimensions where gravity is allowed to propagate [42, 43, 44] in extra 
dimensions. This is not a total coincidence: in models with two large extra dimensions the 
total emissivity of Kaluza-Klein gravitons also scales as Mg~ 4 . We note that due to the high 
power of temperature in the emission rate (3.10) the supernovae constraints are expected 
to be superior to other astrophysical constraints from energy loss mechanisms derived from 
the considerations of red giants, old neutron stars, etc. (See Ref. [41] for further details). 

Similar considerations can be applied to the bremsstrahlung-like emission process N + 
iV — > iV + iV + + 0, where again pairs of are emitted. This process is important 
because the density of neutrons inside the core is rather large. Here, instead of performing 
a detailed calculation which is perhaps not warranted, we use a simple Weizsacker- Williams- 
type estimate for the energy loss. Specifically, we estimate the probability of energy loss in 
a collision of two nucleons to be the product of elastic nucleon-nucleon cross section g^n 
and the probability of an emission of certain amount of energy into 0-quanta by an initial 
or final state nucleon, 

1 T 3 m% f T 



(3 - i3) 
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which leads to an energy loss of 



,2 



1 NN—*NN cfrcfr ~ &NN X X27T 4 M 4 ' ^• i4 i 

where and is the number density of neutrons ~ p/m^. For the relevant range of energies 
<jnn can be taken on the order of 25 mbn [44]. Using the same parameters as before (3.11), 
we arrive at the constraint on M*, 

M* > 15 TeV, (3.15) 

which is very similar to (3.12). With these constraints, we conclude that the effective range 
of (p-force in terrestrial environment (2.9) can indeed be as short as a millimeter. 

Gravitational force constraints. Gravitational force constraints on the model are by far 
the most complicated as they depend very sensitively on the effective distance range this 
force is probed at. 

Could the laboratory measurements of the gravitational force improve over the astro- 
physical bounds (3.12)? First, we look at the exchange by two quanta of (p that does not 
depend on m i n . Such an exchange leads to a 1/r 3 potential, 

v -hw r < 3 - 16 > 

which is limited by recent searches for deviations from the gravitational 1/r behavior at 
short distances [45]. Specifying the constraints on phenomenological coefficient /3 3 from 
Ref. [46] to our model prediction (3.16), we arrive at 

V = .^^i Imm^ with ^ < 13 x 1Q -4 ^ M* > 2 TeV, (3.17) 

which is very close to the astrophysical bounds (3.12) and (3.15). We note again that the 
transition from a 1/r to a 1/r 3 potential that may occur in our model at short distances is 
very similar to the transition expected in theories with two large extra dimensions. 

Constraints from the Yukawa part of (2.12) are somewhat less straightforward to im- 
plement. For a range of O(10~ 2 — 1 cm) in a medium, the constraint on its strength [46] 
specialized to our case with the use of (2.12) takes the following form 

S# S few* 10-, (3,18) 



or 



6L /TeVV / A 2 V i2 

^X(^) X(^r) < 10 12 - 10 
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where we took p ~ 10 g/cm 3 for the density of molybdenum used in experiments of Ref. 
[45]. Perhaps an even more convenient from of the same constraint arises when we trade 
A^ 2 /M* for l/X vac , the range of 0-force in the vacuum, 

Constraints (3.19) and (3.20) do not look intimidatingly stringent, and indeed can be sat- 
isfied by an appropriate choice of M#, 4> m and A 2 . It is also important to note that these 
constraints can be satisfied with a relatively short-range ^-mediated force in vacuum, that 
could be shorter in range than e.g. typical distances within the solar system. 

Clock comparison constraints and constraints on the variations of couplings are of par- 
ticular interest to us in this paper. The now classic comparison of atomic clocks at an 
altitude of R or ut — 10 4 km with clocks on the ground have produced the limit of 2 x 10~ 4 
on possible deviations from predictions of general relativity [47] (as quoted in the review 
[48]), 

< 2 x 10~ 4 x |$(i? ffi ) - <S>{R orbit )\ ~ 5 x 10~ 13 , (3.21) 

where Auu is the extra frequency shift of the Hydrogen maser added to the shift predicted 
by general relativity, and $(r) is the gravitational potential at distance r from Earth's 
center. In our model, the difference between clock frequency on the ground and in orbit 
would receive an additional correction from the difference of coupling constants and masses 
caused by A0: 

uj u a 4 mig p m p L \ 2 J 

where for simplicity we assumed the same scaling for Aqcd an d quark masses with which 
keeps the proton ^-factor g p fixed as a function of distance. The density of a medium 
surrounding the satellite is certainly very low, and thus it is tempting to take A0 2 = <\) 2 m . 
Note however, that if the in-medium range of the force is much shorter than a typical scale 
of a satellite, L sat ~ 1 — 10m, with its average density being M sat /L 3 sat ~ 100kg/m 3 ~ 
O.lg/cm 3 , the scalar field inside a satellite will roll back to its in-medium value, resulting 
in an exponential suppression of A<f>. Therefore, we expect the clock comparison constraint 
to be at the level of 

0^exp(-2L sat /A cff ) x [1 - 2£ e + 4£ F | < 10~ 12 , (3.23) 

noting that the precise amount of exponential suppression would also depend on the po- 
sition of clocks inside the satellite. Although very powerful for M* > 10 8 GeV when the 
exponential factor is of order one, eq. (3.23) is not particularly constraining for an inter- 
esting range of lTeV< M* < lOOTeV, as the value of 0-field inside the satellite would be 
very close to that on the surface of the Earth. 
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Straightforward constraints on 0^ can be deduced from the comparison of coupling 
constants measured in cosmological settings and in the laboratory. Interpreting the null 
results of [17], |Aa/a| < 10~ 5 , in terms of parameters of our model, we obtain the constraint 

|6#m < 2 x lO" 5 . (3.24) 

However, should one accept the criticism expressed in Ref. [18], and interpret the result in 
[4] as a non-zero value of Aa, then one predicts ~ 10~ 5 , and the sign of £p comes 

out to be positive (for positive £p, the low density environments have a larger coefficient in 
front of F* v and therefore a smaller value of a in agreement with [4]). In a similar fashion, 
the indication of a non-zero A(m e /m p ) [19] can be interpreted as a nonzero value for the 
(£e — £p)0m combination. Another interesting possibility occurs for a choice of parameters 
when p ffi ^> A 2 > p qso , where p qso is the average density in a quasar absorption system. In 
this case, not only does one expect a variation in the couplings measured in these systems, 
but one would expect a shift in a from one absorption system to another as function of their 
density 1 . An additional analysis of data in Refs. [4] and [17] searching for a(p) correlation 
might be warranted. 

It is also very important to stress that the Oklo constraint on Aa does not carry any 
weight in our model. Indeed, the Oklo phenomenon obviously occurred in large density en- 
vironment, which means that = m with good accuracy back at the time the Oklo reactor 
was active, as well as it is now. As to the constraints from meteorites [22], similarly to eq. 
(3.23), one expects an exponential suppression of the effect by the density of meteorite, as 
long its size is larger than the 0-field penetration length, L met > X e ff- 

Cosmological constraints. Cosmology can constrain the presence of new degrees of free- 
dom in the Universe. Big Bang Nucleosynthesis can in principle impose a constraint on a 
number of new relativistic degrees of freedom that carries a comparable amount of entropy 
as photons or neutrinos. It is very easy to see, however, that even if is initially thermally 
excited, its decoupling occurs well before the neutrino decoupling because M~ 2 C Gp. 
Since traditionally the BBN constraints are expressed in terms of the number of "new neu- 
trino species", we can immediately conclude that contributes to this number as 4/7 or 
less and thus cannot be ruled out on the grounds of light element abundances [26]. It is in- 
teresting to note that the position of during the time of BBN is close to m , as the energy 
density of nonrelativistic matter during BBN is comparable to the terrestrial p making the 
BBN sensitivity to Aa far less than (3.24). 

The late time evolution of can be rather uninteresting. The scalar field remains at 
= m until the moment when p drops below A 2 , after which it settles towards = 0. 

1 In this sense, the observed scatter in a may be real and due to a local environmental effect. In 
this case, a variation of the density within the absorber. A related solution (which does not require new 
physics) accounted for the scatter through another environmental effect, namely the variation of the isotopic 
abundance of Mg [49] . 
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Since M* <C M P i, the vacuum mass m = A 2 /M* is much larger than the Hubble scale at 
which p ~ A 2 and any oscillations around the minimum are very efficiently damped. Thus, 
the transition from = m to to = was completed when p ~ A 2 ~ eV 4 at redshifts 
2 ^ io 3 — 10 4 . The amount of energy released in such a transition is ~ 0mA 2 , and has a 
negligible effect on the expansion history, as 0^ is constrained to be much less than 1. 

An interesting possibility emerges when A 2 is comparable to the matter energy density at 
redshifts z ~ 1. Then, the evolution of the coupling constants in time towards the minimum 
at = would occur on cosmological scales, while remaining frozen at ~ m on galactic 
scales where the matter energy density is larger than the average cosmological energy density 
by many orders of magnitude. This choice of parameters serves as a concrete realization of 
scenario proposed in [50], where it was suggested that the cosmological evolution of a may 
proceed in some sense independently from the evolution of a in higher density environments 
(see also [51] for a related discussion of AGn as a function of cosmological environment). 
In our model, with A 2 ~ p c , cosmological evolution of would not proceed indefinitely into 
the future but only up to the moment when reaches its low-density minimum at zero. 
Once again, this behavior is possible only because of the increased coupling strength to 
matter, M* < M Ph 



4 Nonlinear models 

In this section we re- introduce nonlinear corrections to 14// (0). In fact, several of such 
models were discussed in the "chameleon" literature [35]- [38], where potentials V e // = 
V(<p) + 0p with different choices of V(4>) were extensively investigated. As we saw earlier, 
the quadratic coupling of to matter allows one to escape many stringent constraints, and 
we intend to keep this feature in this section. To make our discussion more concrete, we 
limit the form of the potential V(<p) to have only even powers of and require the coupling 
to matter to have the same extremum, 

V eff = ^A 4 4 + i(pA 2 + gp)0 2 + A . (4.25) 

Such a Lagrangian may result, for example, from a discrete symmetry — > —0. To ensure 
overall stability, we must choose A 4 to be positive, p and q are taken to be ±1, and we 
maintain A 2 > 0. Another way to interpret the model (4.25) is to say that we allow for 
"tachyonic" value of the mass of in (2.8) but ensure an overall stability of the potential 
by introducing positive 4 contribution. The cosmological consequences of late time phase 
transitions were explored in [52]. 

The most interesting model of this type arises from the choice of parameters that allow 
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for the spontaneous breaking of the discrete symmetry, p = —q. 



P=-l,? = +l =► ( (0 ?-i7a ** '>* = A * . (4.26) 

' = A 4 *(A 2 - p) for p < p c = A 2 . 



It is of course tempting to choose the critical value of matter density p c in excess of the 
average cosmological energy density, so that on average there is a breaking of the discrete 
symmetry in cosmological environments. However, in matter overdensities of sufficient 
spatial extent the symmetry is restored, leading to the erasure of the cosmological vacuum 
expectation value ((f)). The difference between masses and couplings between dense and 
rarefied environments are given by 

Or-Od p-A 2 m r -m d p - A 2 

In the broken phase (low-density environments), there is a linear coupling of the Higgs 
field 0— (<p) to matter, that creates a contribution to the Newtonian force between two test 
particles, 

m 2 ( 9M 2 A — n\ 

U(r) = G N -^ 1 1 + exp(-m cg r) X ^ 2 ^ P j for p < p c = A 2 , ; (4.28) 

"4 = ±=£. (4.29) 

Near the "phase transition", p ~ A 2 , the range of the force mediated by <p becomes infinite. 
In the unbroken phase, only the quadratic coupling of <p to matter survives, which as we saw 
before, significantly relaxes all constraints on the parameters. In other words, the constraint 
(3.17) is still operative while (3.19) is no longer applicable, provided that A 2 < p e . 

As a side remark, we note that there is another interesting spin-off of (4.25) with the 
following choice of parameters: p = — q = 1 and p star < A 2 < p core , where p core is a typical 
density inside a stellar core during a supernovae explosion. In this set-up the scalar field 
keeps its zero expectation value everywhere except for the extremely dense environments 
where it is allowed to roll to ((f) 2 ) ^ 0. Therefore, for this choice of model parameters, there 
are no consequences for the terrestrial and cosmological tests of Act. There is, however, 
an interesting possibility that the nucleon mass experiences a shift during the supernova 
explosion, which in turn may influence the energetics of the explosion and affect the total 
luminosity (see e.g. [53] where the linear change in strong coupling is discussed in connection 
with supernova explosions). If the couplings of neutrons and protons to 2 are different, 
this model can enhance the environmental dependence of the total SN type la luminosity, 
thus affecting the accuracy with which cosmological parameters can be extracted from 
supernovae data. The detailed discussion of such possibilities falls outside the scope of our 
paper. 
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5 Discussion 



On a possibility for new tests of a(p) and m{p) 

One possibility to search for the environmental change of masses and couplings caused 
by a change in density is to try and recreate a low-density environment in the laboratory. 
The best quality vacuums available today achieve a density at the level of 10 5 particles/cm 3 , 
which creates a matter density comparable to eV 4 . Taking p 2 ^> A 2 ^> p±, and using the 
generic solution (3.2), we calculate the resulting shift of and the change in the coupling 
constant between the center and the walls of a spherical chamber of radius R, 



a(r = R) — a(r = 0) 



a 



^G6) 4 fOT R/Koe^l, 

(5.30) 

1 for R/\ vac >l. 



Similar changes will be experienced by masses of particles. Notice that the parametric 
dependence of (5.30) is very similar to (3.20), and plugging A 2 and <p m that saturate this 
constraint we find that 



q(r = R)- a(r = 0) ^ « ! 
a \ M* J 



(10- 17 - 10- 16 ) x (^) 4 for R/X vac « 1, 

(5.31) 

10-15 _ io-i4) x (W) 4 for R/Xvac > 1. 

These shifts are extremely small, but perhaps are not so far away from the modern ca- 
pabilities of frequency measurements that can be sensitive to the relative shifts as low as 
10~ 15 [54]. Further gain in sensitivity can be achieved with substances that have a large 
enhancement factors that connect Au/u with Aa/a [55]. 

To put (5.31) in perspective, we compare this with the result obtained for Aa/a between 
two points separated vertically by ~ lm in the original Bekenstein model with a massless 
scalar field. In that model, the coupling of <fi to matter is linear and very weak, and the 
equivalence principle tests require that (p(Mp\/M*) 2 < 10 -6 (see e.g. Fig 2 of [8]). In 
the Bekenstein model, Aa/a = CfA0, where A(f> is the difference in <fi over lm above the 
surface of the Earth. The solution for is given by 

* = W' (5 ' 32) 

where ( m is the coupling of the scalar to matter (this can be obtained easily from the 
massless limit of the solution given in [39]). In the Bekenstein model, this is absent in the 
Lagrangian but is induced by the coupling of the scalar to the electromagnetic field in the 
nucleon. Roughly, one finds ( rn ~ 10~ 4 Cf [8]. Therefore A<p is simply 

V 3Mf 3 \MJ K ' 
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where h = lm is the vertical separation of the two points. Putting all factors together one 
arrives at 

— (lm) ~ (g x lm) x 1(T 4 C| (^r) £ 10_26 Bekenstein model, (5.34) 
a \M»/ 

where g is again the free-fall acceleration at the Earth's surface. This is well below any 
detection sensitivity for the foreseeable future. The difference between (5.31) and (5.34) is 
an enormous factor of nine orders of magnitude that can be traced back to the fact that 
equivalence principle is checked at macroscopic distances far better than gravity at distances 
under 1 mm. Chameleon models do not allow one to make a simple universal estimate of 
an allowed shift because they depend in a crucial way on the form of the self-interaction 
potential. Widely discussed chameleon models would not allow for a large spatial variation, 
even if one departs from the unviersal coupling of chameleon field to matter. This is again 
related to the extremely tiny values of chameleon field relative to the Planck mass, and 
0(1/Mpi) size of the chameleon coupling to matter. 

Estimates for A(m e /m p ) similar to (5.31) involve a different factor, £ e — ( p = £ e — 1. It 
turns out that this factor can much larger than unity by up to two orders of magnitude. 
The difference comes about due to relaxed astrophysical constraints on £ e /M 2 . Indeed, the 
rate of e + + e~ — > 00 is additionally suppressed relative to (3.10) by a factor of (m e /T) 2 , 
thus allowing for a much larger value of £ e . 

If a nonlinear model of the scalar field such as that described by Eq. (4.25) is realized, 
an artificially created underdensity may lead to the shift of couplings according to (4.27). 
For the shift to occur, one needs p r<R < A 2 < p r >n and R > M*A 2 1 ^ 2 , which would allow 
for the broken phase to be created within r < R volume. 

Astrophysical checks. The possibility of an environmental dependence of coupling con- 
stants calls for new tests of Aa(p) and A(m e /m p )(p). In Bekenstein-type models where 
couplings evolve in time, there is a clear gain from testing a at maximally available red- 
shifts. On the contrary, in the models discussed here, there is no gain in the size of the 
effect at large redshift for tests of Aa(p), and therefore one could search for new test sites 
within our galaxy. Atomic and molecular absorption lines in the interstellar medium could 
be an example where extremely narrow lines can be detected. Another beneficial aspects 
of testing for variations within our galaxy is the possible access to heavier elements where 
relativistic effects are significantly enhanced. It remains to be seen whether O(10~ 5 ) sen- 
sitivity for Aa/a achieved in the QSO absorption spectra can be improved upon using 
galactic lines. 
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6 Conclusions 



We have demonstrated that a scalar coupling to matter can be much stronger than the 
gravitational coupling in Damour-Polyakov type models. The quadratic nature of cou- 
pling to matter allows one to escape the most prohibitive astrophysical and gravitational 
constraints, as only pair-production or pair-exchange of 0-quanta are allowed. The envi- 
ronmental dependence of masses and coupling constant can come about from the shift in 
the expectation value of between dense and rarefied environments. We have shown that 
such shifts could be at a detectable level, and could indeed be probed with astrophysical 
tests within our galaxy and in laboratory clock comparison experiments. 
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